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O ' Abstract 



We study a phenomenological ansatz for merging next-to-next-to-leading order (NNLO) cal- 
>^ . culations with Monte Carlo event generators. We reweight them to match bin-integrated NNLO 

^P \ differential distributions. To test this procedure, we study the Higgs boson production cross-section 

at the LHC, for which a fully differential partonic NNLO calculation is available. We normalize 
PYTHIA and MC@NLO Monte Carlo events for Higgs production in the gluon fusion channel to 
reproduce the bin integrated NNLO double differential distribution in the transverse momentum 
^ . and rapidity of the Higgs boson. These events are used to compute differential distributions for the 

photons in the pp — > if ^ 77 decay channel, and are compared to predictions from fixed-order per- 
turbation theory at NNLO. We find agreement between the reweighted generators and the NNLO 



}-{ ' result in kinematic regions where we expect a good description using fixed-order perturbation the- 

ory. Kinematic boundaries where resummation is required are also modeled correctly using this 
procedure. We then use these events to compute distributions in the pp ^ H ^> W'^W'^ -^ l^l^vv 
channel, for which an accurate description is needed for measurements at the LHC. We find that 
the final state lepton distributions obtained from PYTHIA are not significantly changed by the 
reweighting procedure. 



I. INTRODUCTION 

The search for the Higgs boson is a main objective of the LHC physics program. The 
ATLAS and CMS detectors are designed to detect a Higgs boson in the mass range from 
about 100 GeV up to at least 600 GeV. During the last 15 years, many Higgs boson signatures 
have been studied. For a detailed description of this effort we refer the reader to Refs. [l|, 
y, y, y . If ATLAS and CMS function as designed, a discovery of a Standard Model Higgs 
boson over the entire mass range can be expected with luminosities of about 30 fb~^. 

It is interesting to study how well the mass, width, and couplings of a particle with the 
properties of a Higgs boson can be measured at the LHC, and how well these measurements 
can discriminate between the Standard Model and its viable extensions. The potential 
statistical accuracy for such measurements is usually assessed by computing experimental 
efficiencies using leading order (LO) parton shower Monte Carlo generators. However, these 
efficiency estimates cannot be expected to be accurate, especially in complicated signatures 
such a.s pp ^ H ^ W^W~. 

It has been found that corrections beyond LO are particularly significant for Higgs boson 
production in the channel gg — >■ H. The next-to-leading-order (NLO) QCD corrections [5,^ 
increase the cross-section by a factor of about 1.7-2 as compared to the LO result. A few 
years ago, the inclusive cross-section was computed in the large rritop limit with NNLO 
accuracy [7|,|8|,|9|. The new corrections increase the NLO cross-section by another 20 — 25%. 
Very recently, some threshold-enhanced N^LO terms were also computed 10|, changing 



the NNLO result by less than 5%. These computations show that the Higgs boson cross- 
section can be reliably estimated only after many orders in perturbation theory have been 
considered. It is therefore necessary to account for higher order corrections in a realistic 
analysis of the Higgs boson signal. 

Recentl y, a program FEHIP that describes Higgs boson production in gluon fusion was 
developed julll^l- FEHIP computes the cross-section and fully differential distributions for 



Higgs boson production at NNLO in QCD. Arbitrary cuts can be imposed on partonic jets 
and on the decay products of the Higgs boson. FEHIP does not have an implementation of a 
parton shower and a hadronization algorithm. This creates a few shortcomings, since it is not 
possible to apply cuts at the hadron level or to generate events for a detector simulation. In 
addition, regions of phase-space close to kinematic boundaries can not be described reliably 
in fixed-order calculations. This feature manifests itself in large perturbative corrections at 
special kinematic regions, such as the low Higgs p± region. 

To overcome these characteristic problems of fixed order perturbative computations, the 
resummation of soft gluon effects to all orders in perturbation theory must be performed. 
This resummation may be obtained using analytical techniques |l3L Il4| . Alternativel y, th ese 
effects are also included in parton shower Monte Carlo programs. Novel approaches [13, UM, 
Il9l . |20| merge cross-sections computed in fixed order perturbation theory with these LO 
event generators, such as PYTHIA [15[ and HERWIG [16|. Very significant progress has 
been achieved, and the pioneering Monte-Carlo event generator MC@NLO [131 combines 
consistently NLO perturbative calculations with HERWIG for a number of processes at 
hadron colliders. Unfortunately, no method exists which merges parton shower algorithms 
with NNLO partonic cross-sections consistently. This is desirable for processes with large 
perturbative corrections. 

It is possible to incorporate NNLO corrections into realistic analyses of experimental 
signatures in an approximate way by multiplying probabilities of events in a parton shower 



Monte Carlo simulation by so-called ii'-factors. These factors force the Monte Carlo output 
to agree with certain observables computed in perturbative QCD. This technique is called 
event reweighting. The simplest version of this technique is a multiplication of the Monte 
Carlo output by a constant factor so that the total cross-sections computed perturbatively 
and with the reweighted Monte-Carlo simulation agree. 

Re-scaling Monte Carlo output by a constant factor does not guarantee an agreement 
between perturbative and Monte Carlo results for differential distributions, since perturba- 
tive corrections do depend on kinematic variables and vary across the phase-space. A better 
job may be done if the Monte Carlo output and the perturbative calculation are matched 
at the differential level |21,j. A point-by-point reweighting of the Monte-Carlo throughout 
the available phase-space is not possible, since infrared divergences would produce divergent 
weights. We must instead select a realistic set of observables to match, and then check if 
the reweighted simulation gives a reliable prediction for other observables. 

In this paper we study the reweighting procedure for Higgs boson production in the 
gluon fusion channel at the LHC. We match the Monte-Carlo output of both PYTHIA and 
MC@NLO to distributions that depend only on the Higgs boson kinematics, which is a simple 
and obvious way of reweighting the Monte Carlo output. We match to a double differential 
distribution in the Higgs boson transverse momentum and rapidity. This distribution is 
chosen both for its simplicity and because it allows us to decouple the Higgs boson decay 
chain from the Higgs boson production. However, the kinematics of accompanying QCD 
radiation is totally ignored in the reweighting process. This ignorance is not a problem 
if hadronic radiation is treated fairly inclusively by cuts applied to a process of interest; 
however, if a detailed description of the hadronic radiation becomes relevant, the reweighting 
procedure may lead to inaccurate results. A particular example of a situation when this 
happens is the jet veto on transverse momenta of hadronic jets; we discuss it in detail in 
Section im 

The paper is organized as follows. In Section |H] we introduce the reweighting procedure 
and discuss in detail differential distributions in the reaction pp ^ H + X. We first study 
the reweighting procedure at NLO by comparing the fixed order result with PYTHIA and 
MC@NLO. We present an example in which the reweighting procedure fails to produce 
accurate acceptances: when a jet veto is imposed on the transverse momenta of extra QCD 
radiation. We explain how this problem is ameliorated at NNLO. 

We then apply the reweighting approach to estimate the NNLO effects for the channels 
pp — i> if — i> 77 and pp ^ H ^ W^W~ — >■ l~^i>l~u. We first reweight PYTHIA and 
MC@NLO events in the pp ^ H ^ 'j'y channel. We compute the accepted cross-section 
and differential distributions which have a potential discriminating power from the di-photon 
irreducible background [12, |25|, |26[ : the average pm = ip± +P±)/2 transverse momentum 



distribution and the pseudorapidity difference r]* = \ri'^^ — r^^^l /2 of the two photons. We 
find an excellent agreement between the reweighted PYTHIA and MC@NLO events for all 
observables. The di-photon channel is a testing ground for the reweighting procedure, since 
we can compare the results with the NNLO predictions of FEHIP for the same observables. 
We find that accepted cross-sections agree better than 1%. The Pm, f]* distributions also 
agree very well away from kinematic thresholds. Near these boundaries, they reproduce the 
correct resummed behavior of the parton-shower Monte Carlo simulations. 



We next study the pp ^> H ^> W^W~ -^ l^l^vu channel. Ref. [2l| already employed a 
reweighting technique in order to study the effect of perturbative corrections in this channel 
by matching the PYTHIA output to the resummed pj_ spectrum of the Higgs boson |l4| . 



An optimal set of cuts for isolating a Higgs signal in this channel was introduced and 
studied in |22|. A study of this channel with higher-order QCD corrections included was 
presented in J27|, where both the signal and qq — >• WW background were included using a 
reweighting technique. The analysis in Section IH] of the accuracy achievable by reweighting 
PYTHIA for processes with a jet veto cut implies that at least shapes of distributions 
can be predicted reliably. Hence, we calculate various lepton distributions in the reaction 
pp ^ H ^ W^W~ — > l^l'vu. Interestingly, the reweighting turns out to be largely 
irrelevant for these distributions and the prediction of reweighted PYTHIA and standard 
PYTHIA agree very well. 



II. THE REWEIGHTING TECHNIQUE FOR pp -> H + X 

A. The reweighting procedure 

The cross-sections computed with generator G = {PYTHIA, MC@NLO} for the process 
pp ^ H + X are 

a^ = E/rfn^£(fe})o„(fe}), (1) 

where we sum over all final-state multiplicities m, and integrate the events /^ over the phase- 
space variables dUm of alH < m particles in the final state. The function Om selects the 
kinematic configurations to be accepted in the measured cross-section. The events depend 
implicitly on the renormalization and various factorization scales. 

The simplest observable is the total cross-section af^^i, corresponding to Om {{Pi}) = 1- 
It is a well-known fact that standard event generators fail to predict total cross-sections 
reliably. As an example we set the mass of the Higgs boson to itih = 165 GeV, and 
the renormalization and factorization scales to ^j^r = fip = ''^h/2. We use the generators 
PYTHIA version 6.325 with a Q^ ordered parton shower and MC@NLO version 3.2 For 
PYTHIA we use the MRST2001 LO set of parton-distribution functions, while for MCONLO 
we use the corresponding NLO set. In MC@NLO, the scale choice fJ,R = fip = '"^_ff/2 is also 
used, instead of the default setting. The resulting PYTHIA and MC@NLO cross-sections 
are 

^PYTHIA ^ 12.20 pb ai^,^®^^« = 23.92 pb. (2) 

The corresponding fixed-order NNLO cross-section is 

^NNLO ^ 27.78pb. (3) 

The large differences between the PYTHIA, MC@NLO and NNLO cross-sections reflect the 
fact that the NLO and NNLO perturbative corrections are very significant. 

A consistent method for merging fixed-order perturbative calculations and parton-shower 
algorithms is only formulated at NLO in perturbation theory, and is implemented in 
MC@NLO. A similar procedure beyond NLO is not yet available. Nevertheless, we would 
like to incorporate the large perturbative corrections into the event generators. In this pa- 
per, we adopt a pragmatic approach to solve this problem. We multiply the integrand in 
Eq.(^ with a function K'^, 

a^(«) =^|rfH^/^(fe})ir^(fe})0„(fe}), (4) 



in order to reweight the events f^, 

fG _^ fRiG) ^ fGr^G 
J m J m J m 

The reweighting factors K^ model the effect of higher order corrections through a certain 
order in perturbation theory. We determine the factors K*^ by requiring that Eq.(j3]) repro- 
duces the fixed-order perturbative results for selected distributions, 

(y (Cspccial) = O-^^ (Cspecial) • (5) 

We emphasize that Eq. (Q is an approximate ansatz to describe effects of higher order cor- 
rections in the absence of a rigorous treatment. Strictly speaking, higher order corrections 
do depend on parton multiplicities. For example, pp ^>- H + jets is renormalized differently 
compared to, pp —>■ H + 1 jet. This feature is ignored in Eq.(@]), where the reweighting factors 
Kq do not depend on the multiplicities m. A more detailed version of reweighting would 
not be universal, because matrix elements with fixed multiplicities of partons are divergent 
in perturbation theory. Independent renormalization of events with different multiplicities 
has to depend on a globally-defined set of cuts, e.g. the jet finding algorithm. This in- 
validates the unweightedness of events, the single most important feature of parton shower 
Monte Carlo event generators. We will see the errors in the reweighting procedure caused 
by neglecting the dependence on parton multiplicities later in this section, when we compare 
reweighted PYTHIA at NLO with MC@NLO. 

Having pointed out the approximate nature of the reweighting procedure, we discuss a 
choice of a suitable distribution for which the agreement of a Monte Carlo generator and the 
perturbative calculation can be imposed. Since, as we discussed in the previous paragraph, 
the reweighting ansatz is unsuitable for resolving the structure of QCD radiation, we use 
the kinematic variables which describe the Higgs boson. Since, up to an angle in a plane 
transverse to the collision axis, the Higgs boson kinematics is determined by its transverse 
momentum p± and rapidity Y, we normalize the events /^ to the magnitudes and shapes 
of the NNLO bin- integrated double differential distributions in Y and p±. We expect that 
such a normalization renders the events more realistic in predicting other observables of 
the process. Without a technique for combining NNLO results with parton showering in 
the spirit of MC@NLO, this is the best way we have of combining these calculations with 
event generators. Note that we are not changing the properties of the radiation produced 
by the Monte Carlo generators. We are only changing the normalization of these events to 
reproduce certain distributions. The reweighted generators therefore do not better describe 
events with multiple hard radiations. 

We choose 

if P± e [pi,pi+^] and Y e [Y\Y'+'^] 1 
otherwise, J 



'-^special 

and define the /^-factors as 



(6) 



Kt; := K^iiPf}) = -r\ iipi. e kX^' and Y G y\ F^+^ , (7) 



where Actj^- ' are the accepted cross-section computed at fixed order perturbation theory 
and with the generator G, respectively. The values of the bin boundaries p'^ and F* are 
chosen in such a way that they capture the shape of the Higgs p± and rapidity distributions 



and span the allowed kinematic range for Y and p±. In what follows we always set the 
renormalization and factorization scales to fiR = A^f = iTiMf^, since this choice is known to 
yield a perturbative series with faster convergence J8[ . 
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FIG. 1: Distributions of the Higgs transverse momentum (left) and rapidity (right). 



We point out that at NLO the p± and Y variables completely constrain the kinematics 
of the process pp ^ H + X. At NNLO, new tree level processes pp -^ H + i + j with two 
partons in the final state require additional phase-space variables in order to determine the 
kinematics of partonic radiation. Our reweighting functions K*^ are independent of these 
additional variables. This may create a problem in case there are important cuts on the 
hadronic radiation. 

The choice of the bin boundaries in p^ of the Higgs boson is a subtle issue. Some of 
the standard cuts that we need to apply in Higgs boson production processes have a strong 
dependence on the px_ of the Higgs boson. For example, in the WW decay channel it is 
important to have a phenomenologically realistic model for the low and intermediate p± 
region since after cuts most of the signal comes from the region of low Higgs p±. The 
Higgs p± distribution in this region is not correctly described in fixed-order calculations. 
Logarithms of the form logp±/mH become large and require a resummation. Nevertheless, 
fixed order calculations for cross-sections integrated over pj_ of the Higgs boson are still 
viable, provided that the integration region is sufficiently broad. 

In Fig.^we show the pj_ distributions for the fixed-order NNLO calculation, PYTHIA and 
MC@NLO. We observe that the perturbative NNLO result breaks down at small p±. The 
p± spectrum of PYTHIA is peaked at lower p± than MC@NLO. The most reliable spectrum 
at low p± is obtained with resummation |lj|. To avoid problems associated with the low-p_L 
region in fixed order perturbative calculations, we choose the first p± bin, [p*} = 0,p^], 
in Eq. © to be sufficiently broad by taking p]_ = 25 GeV. Therefore, for p± < 25 GeV, 
we reweight all events with a uniform factor, maintaining the shape of the p± distribution 
provided by the generator G. Above 25 GeV, we trust the shape of the perturbative result 



and reweight in bins of 5 GeV 



p'[ = 0, pi = 25 GeV, pi = (25 



1)5) GeV, 



and 

y^ = 0.5(j-l),j = 1...9. 

Note that this reweighting procedure leads to a discontinuity at p± = p]_ in the reweighted 
p± spectrum computed with the generator G. The choice of p]_ is ambiguous; however, it 
turns out that this ambiguity is largely irrelevant in practice. In what follows we take the 
first bin in p±_ to be [0 — 25 GeV], unless explicitly stated otherwise. 

At this point it is worth investigating if just a single, constant i^-factor is sufficient for 
accurate reweighting. To do so, we investigate the dependence of the reweighting factors 
K'^ in Eq.(j3]) on the p± and rapidity and find that K'^{p±,Y) can vary significantly in 
different rapidity and p^-bins. For PYTHIA, we find K-factors ranging from 1.8 to 3.5, 
while for MC@NLO the /^-factors can vary from 0.7 to 1.6 in bins with a significant number 
of events. For illustration, in Fig. |21we show the reweighting factors for the p± distribution, 
after we integrate over rapidity. We also show the reweighting factors as a function of Y, 
after integrating over p±. The shape of the X-factors in the two variables is not uniform, 
indicating that a naive multiplication with a uniform i^-factor from the total cross-section 
may not be adequate. Having discussed the reweighting technique in general, we now study 
it in the process pp ^ H + X at both NLO and NNLO. 
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FIG. 2: The reweighting factors integrated over rapidity for PYTHIA and MC@NLO as functions 
of the Higgs p± (left). The reweighting factors, integrated over p±, as functions of rapidity (right). 
The inclusive iT-factor for the total cross-section is also shown on both plots as a horizontal line. 



B. The reweighting at NLO and NNLO 



We now apply the reweighting procedure to PYTHIA and MC@NLO to study pp — *■ 
H + X. MC@NLO is a Monte Carlo event generator that accommodates NLO perturbative 





calculations for a wide range of processes. Its important feature is that the parton shower 
is combined with NLO matrix elements consistently, without double counting. PYTHIA is 
an event generator based on leading order matrix elements, so that hadronic radiation is 
generated primarily through the parton shower. 

We begin our study by checking how well reweighting works at NLO. We check how 
well the procedure describes Higgs boson production when a jet veto is imposed. This 
tests whether the neglect of the additional hadronic radiation in our reweighting ansatz is 
problematic. Phenomenologically, this cut is needed in the pp ^ H —> W~^W~ — > l~l'^vu 
channel to isolate the signal from background. The inclusive cross sections for both PYTHIA 
and MC@NLO are given in Eq.©, while a^^^ = 23.99 pb. We impose a jet veto of p_L < 30 
GeV, and we define jets with a cone algorithm using a cone size R = 0.4. The cross sections 
after the jet veto has been imposed are 

PYTHIA; 

. , . _ . ...... MCONLO; ,_. 

0"acclPDj - "^ -iAASi RNLOmyTHIA); ^^ 

NLO. 

The acceptances, defined as the ratios of the accepted cross sections over the inclusive cross 
sections, are 

PYTHIA; 
_ . ...., MC@NLO; , . 

^-^nsn RNLO(PYTHIA); ^^^ 

NLO. 

We observe a very large disagreement, of order 30%, between the acceptances obtained 
using the generators and the fixed order NLO result. What is occurring here is that this 
observable is very sensitive to the properties of the QCD radiation. Multiple partonic emis- 
sions are required to generate the correct jet p± spectrum, and the NLO result contains only 
a single partonic emission. The p± spectrum of this additional parton is generated for the 
first time at NLO, and is therefore not accurately predicted at this order in the perturba- 
tive expansion. We note that reweighting PYTHIA to the NLO result spoils the agreement 
between its acceptance and that computed with MC@NLO. 

To check that multiple emissions are indeed important, we present below the jet mul- 
tiplicites for both PYTHIA and MC@NLO before a jet veto is imposed. We study the 
cross-section both inclusively and with the restriction p^ > 30 GeV, to show that multiple 
emissions are required to obtain correctly even the high pj_ events. We present the fraction 
of events with A^ = 0, 1, 2, 3, 4, or more jets in Table |l] Note that we require a jet to have 
p± > 20 GeV, so events without jets are possible. Over half of the events in the high p^ 
tail coming from PYTHIA and MC@NLO contain multiple emissions, indicating that the 
description of the hadronic radiation coming from the single emission at NLO is unlikely 
to be very accurate. We can also see this by studying the Higgs pj_ spectrum, shown in 
the left panel of Fig. |21 The single hard partonic emission is equivalent in the fixed order 
NLO result and in MC@NLO. The mismatch between them in the high p^ tail is caused 
by showering. The importance of the multiple emissions is made explicit in the right panel 
of Fig. 121 where the p^ from MC@NLO when only a single jet is observed is compared to 
the NLO calculation. The distributions agree very well in the high p^ region for this single 
emission case, again indicating the need for multiple emissions to correctly generate this 
spectrum. 
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Inclusive 


pf >30 




PYTHIA 


MC@NLO 


PYTHIA 


MC@NLO 


A^ = 


0.365 


0.39 


0.055 


0.090 


N = 1 


0.335 


0.345 


0.40 


0.465 


N = 2 


0.18 


0.17 


0.31 


0.275 


N = 3 


0.080 


0.060 


0.15 


0.105 


N = A 


0.030 


0.020 


0.055 


0.040 


N>A 


0.010 


0.015 


0.030 


0.025 



TABLE I: Fraction of events with A^ = 0, 1, 2, 3, 4 or more jets for inclusive Higgs boson production 
and Higgs boson production with p^ > 30 GeV in PYTHIA and MC@NLO. 
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FIG. 3: The Higgs boson p± spectrum for NLO, MC@NLO, and PYTHIA reweighted (left panel); 
comparison of the Higgs boson p± at NLO, and with MC@NLO when only a single jet is observed 
(right panel). 



We next study what happens when we perform the reweighting at NNLO. We use FEHIP 
to obtain these results. The inclusive NNLO cross section is given in Eq.Q. We include the 
reweighting of MC@NLO to the NNLO double differential distribution. The accepted cross 
sections for NNLO and the reweighted event generators are 



cr„ 



:(pb) 




j^NNLO (PYTHIA); 

j^NNLO(MC@NL0); 
NNLO. 



(10) 



The acceptances are 



A 




RNNLO(PYTHIA); 

RNNlo(MC@NLO); 
NNLO. 



(11) 



We observe a much better agreement with the NNLO reweighting. R(PYTHIA), 
R(MC@NLO) and the fixed order NNLO result all agree with the PYTHIA and MCONLO 
acceptances within 6%. The NNLO result contains two partons in the final state, which 
gives a more realistic accounting of the QCD radiation. It also contains the first radia- 
tive correction to the single parton p±_ spectrum. The p± spectrum obtained at NNLO is 
in better agreement with MC@NLO, as seen in Fig. ^ A comparison of the pj_ spectrum 
from the reweighted generators with the resummed p± distribution of J28| is presented in 
Fig. 0] There is good agreement between R^^'"°(MC@NLO) and the resummed calcula- 
tion. R^'^^*-' (PYTHIA) agrees with the intermediate and large p± portion of the resummed 
distribution, while there is a slight discontinuity induced by the first bin reweighting in 
the low p± region. We conclude that even in the presence of significant cuts on the jets 
in the final-state, the simple reweighting of the Higgs boson double differential distribution 
at NNLO describes the acceptances well. In addition, since the NNLO result produces the 
correct normalization and contains drastically reduced scale dependences, we believe that 
reweighting MC@NLO with the fully differential NNLO result of FEHIP provides a very 
accurate prediction for the Higgs boson signal at the LHC 

As a final check, we compute the rapidity distribution of the Higgs boson using FEHIP 
and the reweighted event generators. The result is shown in fig. El We observe that imposing 
the jet veto maintains the matching of this distribution. 
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FIG. 4: The Higgs boson p± spectrum from a resummed calculation, MC@NLO reweighted and 
PYTHIA reweighted. 

Motivated by the success of the NNLO reweighting procedure, we now allow the Higgs 
to decay and study predictions for the pp— >-if + X— >-77-|-X and pp — >• if — > W^W~ -^ 
l~l^vv channels. 
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FIG. 5: The Higgs boson rapidity distribution with a jet veto of 30 GeV apphed for NNLO, 
MC@NLO reweighted and PYTHIA reweighted. 

III. THE DI-PHOTON CHANNEL 

We first use the reweighting procedure to compute the cross-sections for observables in 
the pp—*H + X— >■'-/'-/ + X channel. The standard cuts on the two photons are 

• pf > 40 GeV and p]^ > 25 GeV; 



\r]' 



71,2 



< 2.5 ; 



• -Ehadr < 15 GeV in cones of size R = 0.4 around each photon. 

The two-photon channel is useful for additional checks of the reweighting approach. Because 
the if — > 77 decay is included in the FEHIP program, we can compare observables computed 
with R(PYTHIA) and R(MC@NLO) directly with NNLO resuhs. 

We first compare the various results for the accepted cross-sections after applying the 
standard cuts. We choose a Higgs mass rriH = 120 GeV and set the renormalization and 
factorization scales to Hr = fJ^F = tt^h/'^- For the if ^ 77 branching ratio we assume the 
value Br(ii -^ 77) = 0.002205, which we obtain from HDECAY 24]. We find the following 
cross-sections for the di-photon signal: 





PYTHIA 


MC@NLO 


i?(PYTHIA) 


i?(MC@NLO) 


NNLO 


C^acc [fb] 


36.8 


60.3 


65.3 


66.9 


66.4 



The PYTHIA, MC@NLO and the NNLO results differ significantly, reflecting again the large 
NLO and NNLO corrections. However, the reweighted cross-sections agree within 2.5%, and 
differ from the NNLO result only by -1.7 % for R(PYTHIA) and +1.0 % for R(MC@NLO). 
The effect of the cuts on the accepted cross-section in the two-photon channel is rather 
insensitive to the choice of generator, and the reweighting procedure reproduces the NNLO 
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results reliably. This is an expected result, since in the di-photon decay the experimental 
cuts do not resolve the structure of the hadronic system that recoils against the Higgs boson. 
We next compare the reweighted results and the NNLO predictions for more compli- 
cated observables in the di-photon channel. In Refs [12, |25|, |2y| the distribution of the 
pseudorapidity difference of the two-photons y* = 1/2 \ri^'^ — ri^^\ was proposed as a discrim- 
inator to separate the signal from the prompt photon background. In Fig. [Ulwe present the 
?/*-distribution for PYTHIA and MC@NLO, as well as for the reweighted generators and 
at NNLO. We observe that the reweighted and the NNLO distributions agree reasonably 
well. At the boundary of the kinematic region y* > 0.96, the NNLO distribution is non- 
monotonic. This kinematic region corresponds to a vanishing leading order cross-section, 
and the perturbative result must be resummed to all orders. PYTHIA and MC@NLO do not 
suffer from this problem since parton showers perform such resummations. The reweighted 
generators maintain the resummed behavior at the kinematic boundary and reproduce the 
fixed order result elsewhere. 
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FIG. 6: The pseudorapidity difference y*-distribution for di-photon events. We compare PYTHIA 
and MC@NLO with the reweighted generators and NNLO. 



We now study the average p± distribution of the two-photons, Pm = 1 /2 ip± +P± ) [1^ I2f 
At leading order in perturbation theory the cross-section is zero for p^ > mH/'2. The 
distribution at higher orders retains a characteristic peak at pm ~ mH/'2,. In Fig. [7| we 
show the distribution for PYTHIA, MC@NLO, and FEHIP, and also after reweighting. We 
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FIG. 7: The average transverse momentum pm-distribution for the two photons. We compare 
PYTHIA and MC@NLO with the reweighted generators and NNLO. 



again find a very good agreement between the R(PYTHIA) and R(MC@NLO) results. The 
NNLO distribution agrees very well away from the peak at Pm ~ rnH/2. As expected, the 
NNLO result at the peak is substantially different because this region cannot be predicted 
accurately in fixed-order perturbation theory, and requires resummation. The reweighted 
generators do a reasonably good job at maintaining the appropriate resummed behavior at 
the peak. 

The ratio of the di-photon cross-sections computed with the reweighted and the leading 
order generators is very similar to the NNLO i^-factor for the total cross-section. The 
reweighting of the event generators with a constant factor could also yield realistic results for 
the di-photon cross-section after applying the standard cuts. It is interesting to investigate 
if a constant i^-factor is also sufficient for reweighting differential distributions. We have 
already seen that this would not be satisfactory for the p± and rapidity distributions for the 
Higgs boson in Fig. |21 To investigate this, we compute the effective i^-factors for each bin 
of the y* and Pm distributions. 



Kcihin) 



Aa^(^)(bin) 
AaG(bin) 



(12) 



13 



2.8 
2.6 
2.4 
2.2 
2.0 
1.8 
1.6 
1.4 
1.2 
1.0 



- 


LHC14TeV 




r 


m„ = 120GeV 

— 1 , . r 




1 






•^inc 



0.0 







- 






1.4 


- 






1.2 


- 




J 


F — — ' — 1 . " ' 




1.0 




^ 




0.8 


— 






0.6 







0.2 



0.4 



0.6 

y* 



0.8 



1.0 



1.2 



FIG. 8: The effective /sT-factors as a function of y* for PYTHIA and MC@NLO. 

In Fig. ISlwe show that the effective i^-factors in each bin of the y* distribution do not vary 
significantly from the inclusive iiT-factor. However, in the average photon p± distribution of 
Fig. the effective i^-factors for high pm bins away from the peak are not uniform, although 
the large kinematic variations occur in bins with few events. 

IV. THE W+W- -^ l+l-uu CHANNEL 



In this Section we study the reaction pp —* H —>■ W^W~ —^ l^l^vv. This signal can 
be distinguished from the main background process pp —>■ W^W~ by applying a jet veto, 
requiring a small opening angle between the two charged leptons in the transverse plane, 
and applying some additional kinematic cuts J22ll23j. 

In the remainder of this Section, we present new results with the R(PYTHIA) generator 
in the H — >■ W^W~ -^ l~^l~i>u decay channel. The observables that we consider probe 
the momenta of the final-state leptons. The W^W~ decay chain is not yet implemented 
in FEHIP, and a comparison with NNLO is not possible. In addition, the Herwig event- 
generator, which is a basic component of MC@NLO, does not have an implementation of 
the same decay with full spin correlations. Therefore, we will not present any leptonic 
observables with R(MC@NLO). 

We first present a study of the ii'-factor for the H -^ W^W~ -^ l^l^vv channel. Using 
R(PYTIIIA), we can obtain a description of the Higgs boson rapidity distribution valid to 
NNLO, and examine the effect of this rapidity dependence. The effective i^-factor integrated 
over the whole region after all H —>■ WW cuts are applied is 2.098. If we reweight PYTHIA 
to only the p± spectrum of the Higgs boson, we find an effective i^-factor of 2.02. The 
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FIG. 9: The effective /sT-factors as a function of the pm distribution for PYTHIA and MC@NLO. 

{p±,Y) dependent effective K-factoT is 10% lower than the fully inclusive K-iactoi of 2.28, 
while the iiT-factor coming from only reweighting to the p \ di stribution is 13% lower than 
the inclusive /^-factor, comparable with the results from J28|. The effect of the rapidity 
dependence is therefore less than 3%. 

We now employ R(PYTHIA) to study distribution shapes in pp — > i7 — >■ W^W~ -^ 
l^l^vv. We normalize each result to the integrated cross-section subject to the appropriate 
cuts. Since the distributions studied do not probe the hadronic radiation, we expect them 
to be very well described by R(PYTHIA). 

In Fig ^1 we plot the minimum and maximum transverse momentum distributions of the 
detected leptons for PYTHIA and R(PYTHIA) events. These distributions are characteristic 
of the Higgs signal and can be used to discriminate from the background. We observe that 
the reweighting does not change the shape of distributions. An application of a constant 
i^-factor would lead to the same results. However, the appropriate i^-factor is the effective 
one of 2.098 discussed at the beginning of this section, which is 10% lower than the fully 
inclusive i^-factor. 



V. CONCLUSIONS 

In this paper, we have presented a phenomenological approach to including NNLO cor- 
rections in event generators such as PYTHIA or MC@NLO. Without an extension of the 
MC@NLO procedure to NNLO, this offers the best way of combining parton showering and 
hadronization with NNLO calculations. We study this procedure for Higgs boson production 
at the LHC, since the fully differential NNLO calculation is available in the program FEHIP. 
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FIG. 10: The minimum (left) and maximum (right) transverse momentum of the two leptons 
computed with PYTHIA and R(PYTHIA). 

We reweight the Monte-Carlo events of PYTHIA and MC@NLO to match the bin-integrated 
NNLO double differential distribution in the Higgs p± and rapidity. We then study how well 
distributions of the Higgs boson decay products are predicted by this reweighting procedure. 
We note that the i^-factors that describe the reweighting of both PYTHIA and MC@NLO 
exhibit significant kinematic dependences, so that the use of a single constant i^-factor may 
not be adequate. 

We first study the reweighting procedure for the process pp -^ H + X, without de- 
cays of the Higgs boson. The i^-factors that describe the reweighting of both PYTHIA 
and MC@NLO depend significantly on the Higgs boson transverse momentum, and non- 
negligibly on its rapidity. We test how well reweighting reproduces a fully consistent merging 
of fixed-order calculations with parton showering by reweighting PYTHIA at NLO and com- 
paring to MC@NLO. We find large discrepancies when a jet veto is imposed. This indicates 
that the single parton emission present in the NLO calculation is insufficient to correctly 
describe cuts where the hadronic structure is probed. The reweighting in the presence of 
a jet veto works much better at NNLO, where two partonic emissions are present in the 
final state. This indicates the importance of extending perturbative calculations to NNLO 
in order to obtain a reasonable description of the additional radiation. 

We then examine the decay channel pp ^ H -^ '-yy with all relevant experimental cuts 
included. We find that both the reweighted PYTHIA and the reweighted MC@NLO match 
very well the accepted cross-section as predicted by FEHIP. We next study distributions 
that have been proposed to discriminate between the Higgs signal and the background. Both 
R(PYTHIA) and R(MC@NLO) describe the kinematic distributions well. They match the 
NNLO fixed-order result away from kinematic features, and exhibit the resummation present 
in the event generators near the kinematic boundaries. 

We proceed to study the decay channel pp ^ H ^ W^W~ -^ l^l^vv. It is important 
to understand distributions in this channel at the LHC, since a direct reconstruction of 
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the Higgs boson mass peak is not possible because of the two neutrinos in the final state. 
We can not yet directly compare lepton distributions with the NNLO result, since FEHIP 
does not yet contain the decay H -^ W^W~ — >• l^l'vp. We study lepton and missing 
energy distributions using R(PYTHIA), assuming that it predicts the distribution shapes 
correctly. Since Herwig does not yet contain spin correlations for this channel we do not 
present results for R(MC@NLO). We find that the reweighting induces very small kinematic 
shifts. We study the effective ii"-factor for this channel after all cuts have been applied, and 
find that it is 10% smaller than the inclusive K-{a.cioi. The effect of the rapidity dependence 
on this effective i^-factor is small, about 3%. 

In summary, in this paper we study for the first time the detection efficiency for the Higgs 
boson at the LHC by reweighting parton shower Monte Carlo output to the fully differen- 
tial Higgs boson cross section at NNLO in QCD. Monte Carlo events from PYTHIA and 
MC@NLO are normalized to the NNLO calculation of the Higgs boson rapidity and trans- 
verse momentum distributions. For Higgs boson events with low p^, a constant i^-factor is 
applied, to maintain the resummed shape present in the Monte Carlo simulations. For the 
H -^ WW channel , where a jet veto is applied, we find a small difference of 3% compared to 
a similar reweighting approach using only the transverse momentum spectrum. We conclude 
that the effect of the NNLO Higgs boson rapidity dependence on LHC observables is now 
accurate to the percent level. The dominant remaining theoretical uncertainties affecting 
the Higgs boson search at the LHC are: (1) the scale uncertainty arising from truncation of 
the perturbative expansion at NNLO; (2) the modeling of the low p±_ Higgs spectrum; (3) 
the theoretical uncertainties for the backgrounds to the Higgs signal. We believe that the 
event reweighting studied here is a useful and accurate way of including higher order QCD 
calculations in Monte-Carlo event generators. It allows us to determine the correct normal- 
ization at NNLO after all experimental cuts are included, incorporates the kinematic shifts 
induced by hard QCD radiation at higher orders, and maintains the resummation present in 
parton-shower Monte-Carlo programs near kinematic boundaries. We believe that reweight- 
ing MC@NLO to match NNLO differential distributions gives a highly accurate description 
of the Higgs signal at the LHC. We look forward to its application to other processes of 
phenomenological interest. 
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